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The collective excitation spectra are studied for a multisubband quasi-one-dimensional electron 
gas on the surface of liquid helium. Different intersubband plasmon modes are identified by calculat- 
ing the spectral weight function of the electron gas within a 12 subband model. Strong intersubband 
coupling and depolarization shifts are found. When the plasmon energy is close to the energy differ- 
ences between two subbands, Landau damping in this finite temperature system leads to plasmon 
gaps at small wavevectors. 



The electron system on the surface of liquid helium 
discovered about 30 years ago has provided an ideal plat- 
form to study many-body effects in low dimensions. [1,2] 
Recently, quasi-one-dimensional (Q1D) electron gas em- 
bedded in such structures has been realized in labora- 
tory by bending the liquid helium surface. [3] In some 
sense, this Q1D system is similar to those in semicon- 
ductor quantum wires, nanotubes and metallic chains 
where many-body effects have been extensively studied. 
[4] However, this new system on the surface of liquid 
helium provides us more degrees of freedom to explore 
many-body effects in a ID framework. It is free from 
impurities being much clearer than other systems. Fur- 
thermore, it is of a wider variable range of electron den- 
sities. Several theoretical and experimental studies have 
been carried out on the electron transport and electron- 
electron interactions in this system. [5,6] The electron 
system on the liquid helium surface is achieved at finite 
temperatures and the electron density is usually much 
smaller than that in semiconductor structures. In most 
cases, it is considered as a non-degenerate electron gas 
obeying the Maxwcll-Boltzmann statistics. [7] On the 
other hand, the ID confinement is weak and the gap 
between the ID subbands is of the same order as the 
thermal energy. As consequence, many subbands are 
usually occupied forming a real multisubband Q1D elec- 
tron system. Very recently, dispersion relation of the col- 
lective excitations (plasmons) was studied theoretically 
within the random-phase approximation (RPA) for the 
Q1D electron gases on the surface of liquid helium. [7,8] 
The RPA is believed to be a more reliable approximation 
to study collective excitations of such a non-degenerate 
Q1D electron gas in the classical regime. In these studies, 



[8] a two-subband model was used and particular atten- 
tion was devoted to the dispersion relations of one in- 
trasubband and one intersubband plasmon modes found 
within this model. Notice that, both the RPA and the 
quasi-crystalline approximation, being valid in opposite 
limits of electron density at fixed temperature, result in 
very similar plasmon dispersions indicating that the RPA 
could be correct over a wide range of electron densities. 

This Rapid Communication is focused on the collec- 
tive excitation spectra of the multisubband Q1D electron 
gas on the surface of liquid helium. Special attention 
is paid to the intersubband coupling and effects of the 
single-particle excitations (SPE), i.e. the Landau damp- 
ing, on the plasmon modes at finite temperature. We 
study these elementary excitations by calculating the full 
spectral weight function of this classical electron system. 
The Landau damping induced phenomena are expected 
to be of a new kind provided the electrons are in a very 
different regime from those embedded in semiconductor 
quantum wires. In a Q1D electron system, the inter- 
subband interactions are much stronger than those in 
higher dimensions. Furthermore, temperature effects are 
important because the energy gap between the ID elec- 
tron subbands is comparable to the thermal energy ksT . 
Many subbands can be occupied even for small densities 
and temperatures. As a consequence, the intersubband 
coupling should be treated properly. In the calculations, 
we take into account as many subbands as possible to 
guarantee the efficiency of the method and to make sure 
the obtained plasmon spectra being independent of the 
number of subbands included. 

We consider the same structure for the Q1D electron 
system as in Ref. [7] where the electron mobility was 
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studied. The ID confinement on the liquid helium is 
determined by both the surface curvature radius in the 
y-direction R = 5 x 10~ 4 cm and the so-called holding 
field applied in the ^-direction E± = 3 x 10 3 V/cm. 
This confinement is then approximated by a parabolic 
potential in the y-direction with a confinement frequency 
loq = -J eE±/mR, where m and e are the mass and the 
charge of the electron, respectively. The energy eigen- 
values of an electron in the system is given by E n {k x ) = 
k x /2m + (n — \)uJq, where k x is the electron wavevec- 
tor in the x-direction and (n — ^)ojq is the electron en- 
ergy levels due to confinement in the y-direction, with 
n = 1,2, being the subband index. Here, we consider 
the electron gas being of zero-thickness in the ^-direction 
because the energy gap between the two lowest levels due 
to confinement in this direction is greater than 30 K, [9] 
whereas the confinement energy loq in the y-direction is 
less than 1 K. Another important parameter is the lo- 
calization length 7/0 = k^ 1 = \JTfZnuJ~o of the electrons 
in the y-dircction (we consider h = 1 throughout this 
paper). 

Our calculations show that there is only one observ- 
able intrasubband plasmon mode but many intersubband 
modes in the present Q1D electron system. All the ob- 
servable intersubband plasmon modes are related to the 
first and second subbands denoted by (l,n) and (2, n), 
respectively, where n = 1,2,3,.... Such intersubband 
modes are important since most of the electrons occupy 
states in the two lowest subbands and intersubband inter- 
actions are strong in a one-dimensional geometry. Fur- 
thermore, we find that intersubband coupling between 
the higher subbands and the two lowest subbands af- 
fects significantly the low-energy intersubband plasmon 
modes. A correct theoretical consideration on such a 
coupling turns out to be essential in obtaining even the 
lowest intersubband plasmon modes. In order to con- 
sider this effect properly, we use in this work a 12 sub- 
band model so that the plasmon spectra shown below do 
not change anymore when more subbands are included. 
Moreover, we show that pronounced Landau damping 
occurs when the plasmon energy approaches the energy 
difference ncj between two subbands leading to gaps in 
the plasmon spectra. 

A traditional way to obtain the plasmon dispersion re- 
lations in Q1D electron gases at zero temperature [10] is 
to find the roots of the equation det \Ke{e a /3(q, = 0, 
where e a p is the multisubband dielectric matrix and the 
indices a = and j3 = (j,f) with the subband in- 

dices i, i',j and j'. However, this method cannot provide 
us complete information of the plasmon excitations in the 
present system. There are two main reasons for that: (i) 
it is difficult to figure out where the plasmon excitations 
are Landau damped by SPEs because Im{e a p(q, u))} ^ 
in the whole ui-q plane at finite temperatures; and (ii) 
the thermal fluctuations in the system might easily pop- 
ulate several subbands even for small densities. There- 



fore, intersubband interactions can be strong so that a 
reasonable amount of subbands should be included in the 
calculation. The equation det |Re{£ Q/ 3}| = yields many 
roots but does not provide information of the relative 
importance of each one. Thus, one cannot distinguish 
which roots correspond to the plasmon modes. 

A more reliable way to study the collective excitations 
in such a multisubband system at finite temperatures is 
to calculate the so-called spectral weight [11] 
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The peaks of this function give information of the plas- 
mon excitations and, from their position, we can obtain 
the dispersion relations of the plasmon modes. By show- 
ing all excitation modes through S(q,w), one provides 
a very efficient guide of what should be observable in 
the experiments. These observations are certainly de- 
pendent on external probes (e.g., light polarization). Eq. 
1 also involves the ID non-interacting irreducible po- 
larizability function U^(q,tu). We remember that the 
Maxwell-Boltzmann distribution function is used in cal- 
culating Hg(q,u)), with 5 being a phenomenological con- 
stant which is responsible for broadening of the energy 
levels E n {k x ) mainly due to ripplon and evaporated he- 
lium atom scattering on the surface. The quantity S(q, u>) 
is directly related to the optical (such as inelastic light 
scattering spectra) and transport (such as conductivity) 
properties. The multisubband dielectric matrix function 

£ a p{q,u>) = 5 af 3 - V a p(q) ■ U 5 p (q,w) 

is written within the RPA, with V a 0(q) being the 
Coulomb electron-electron bare interaction. 

The symmetric confinement potential in the y- 
direction leads to the electron-electron Coulomb inter- 
action V a p (q) = when i + i' + j + j' is an odd num- 
ber. Consequently, the dielectric matrix elements (both 
the real and the imaginary parts) E a p(q,u>) = for 
i + i' + j + j' = odd. The dielectric matrix can then be 
decoupled into two submatrices e^S n (q, oj) and e°^ d (q, w) 
with both i+i' and j+f being even and odd numbers, re- 
spectively. The even (odd) dielectric submatrix involves 
only intersubband electron-electron interaction for one 
electron from subband i to subband i' with i + i' = even 
(i + i' = odd) and the other from subband j to f with 
j + j' = even (j + j' = odd). As a consequence, the 
spectral weight function can be treated separately in two 
parts S(q,u) = S cvca (q, uj) + S odd (q, lo). 

In order to understand the damping induced phenom- 
ena in the present system, we calculate first the spec- 
tral weight of the single-particle excitations S sp (q 1 uj) = 
— J2jj> I m [n|j/ (<?, w)] • It is straightforward to analyt- 
ically obtain the imaginary part of the polarizability 
ILp(q,u)) at finite temperatures. [8] Considering the level 
broadening effects, it is given by 
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Im [II*.] =AT e {exp[-i/T] I s (£<+>)- 

expH/T] (47)) } /qV^T [1 + coth (1/2T)] (2) 

where N e is the total electron density, T the temperature, 
q the wavevector in the ^-direction, and 



-I 



dfc. 



exp[-(fc a; +g) 2 ] 
fc? + A 2 



(3) 



with = [w + (i - j) ± g 2 ] /2<7\/Tand A = S/2qVT. 
The wavevectors fc x and g are written in unit fco , while the 
frequency ui, the broadening constant 5 and the tempera- 
ture T are in unit ojq. Notice that, in the limit S — > 0, the 
integral / 5 (£) = 7rcxp(— £ 2 ), which is a Gaussian func- 
tion of £. We take the broadening parameter S = 10~ 2 wo 
throughout this paper. This value corresponds to realis- 
tic electron mobilities due to ripplon scattering. [7,12] 

In Fig. 1 wc show S sp (q,uj) as a function of u) for 
different values of the wavevector q. Here, the temper- 
ature T = 0.4 K, the electron density N e = 10 4 cm -1 
and the holding field E±_ = 3 x 10 3 V/cm. The en- 
ergy gap between the quantized levels at this holding 
field is uj ~ 0.784 K which leads to electron occupation 
up to the second and the third subbands of about 14% 
and 2%, respectively, in relation to the first one. For 
the sake of consistency, we plot S sp (q,u)) in two parts: 

(9.^) = -Ej+/=even Im [ U jj' h^)] ( Solid Curves) 



ceven 



and S° dd (q,uj) = -£ 



Im [n^., (q, lu)] (dashed 



j+j'= odd ' 

curves). For small wavevector q, the spectral weight of 
the intersubband SPEs are mostly described by Gaus- 
sian functions centered at energies 10 = nwo, where 

n = 1,2,3, On the other hand, the intrasubband 

SPE are represented by the less weighted peak seen in the 
lower energy part of the spectra. This peak is mostly due 
to excitations in the first subband. Fig. 1 show us that 
SPEs 'induced effects should be significant at u = nuj 
for small q. 

We show in Fig. 2 the spectral weight (a) S even (q,u>) 
and (b) S odd (q,ui) for different electron densities N e at 
a small fixed wavevector q = O.lfco- The temperature is 
taken as T = 0.4 K. The density varies from N e = 1 x 10 3 
cm -1 (the lowest curve) until N e = 1.45 x 10 4 cm -1 
(the top curve) with a difference 1.5 x 10 3 cm" 1 . For 
N e = 1 x 10 3 cm _1 , only the first subband has a signif- 
icant electron density so the plasmon modes related to 
this subband can be observed. This help us to identify 
the different peaks corresponding to the intrasubband 
mode (1,1) and intersubband modes (1,3), (1,5), and 
(1,7) in Fig. 2(a). The peaks in Fig. 2(b), at the same 
electron density, are due to the plasmon modes (1,2), 
(1,4), and (1,6). As the density increases, these peaks 
shift to higher energy because the depolarization shift is 
enhanced. Meanwhile, the intersubband modes related 



to the second subband appear as indicated in the fig- 
ure. We also observe that the Landau damping becomes 
stronger when the plasmon mode (l,n) approaches to 
the frequency [n + l)o>o where the single-particle exci- 
tations are of high intensity as shown in Fig. 1. The 
insets show the energy position of the peaks due to the 
intersubband plasmons as a function of the density N e . 
The gaps result from the Landau damping and are clearly 
seen around the energies uj = 4a> an d 6o>o in Fig. 2(a); 
and u) — 3ui n and 5uo in Fig. 2(b). 

In Fig. 3, we plot the spectral weight function for dif- 
ferent q values corresponding to those in Fig. 1 with 
N e = 10 4 cm- 1 and T = 0.4 K. Figs. 3(a) and 3(b) show 
S cvcn (q,oj) and S odd (q,Lj), respectively. The dispersion 
relations of the plasmon modes obtained from the peak 
position are given in the insets. We also see small peaks 
in the full spectra due to single-particle excitations indi- 
cated by the open-dots in the insets. The lowest branch 
(open-dots) represents the intrasubband SPEs. Similarly 
as in Q1D Fermi liquid electron systems at zero temper- 
ature, with increasing q, the spectral weight of the in- 
trasubband (intersubband) plasmon mode increases (de- 
creases). We also observe the Landau damping induced 
gaps appearing around w = nuj at small q where the 
SPEs are strong. Furthermore, the peaks due to the 
plasmon excitations become wider and lower at larger 
q because S sp (q,u)) tends to be a more uniform function 

of CO. 

For the sake of completeness, we analyze in Fig. 4 
the spectral weight S(q,u) as a function of the energy 
uj and for several temperatures. Here, the density N e 
= 10 4 cm -1 and q = O.lfco. From the bottom to the top 
curves, the temperature increases from T = 0.2 K to 2.0 
K with a step of 0.2 K. As the temperature increases, 
the peaks of the intersubband plasmon modes (l,n) re- 
lated to the first subband shift to lower frequency while 
those modes (2, n) related to the second subband shift 
to higher frequency. Such a change is mainly induced 
by the redistribution of the electron density in different 
subbands. As temperature increases, the electron den- 
sity of the first subband decreases but that of the higher 
subband increases. The decrease of the electron density 
in the first subband results in a decrease of the depolar- 
ization shift of the related intersubband plasmon modes 
(l,n). The insets show the energy position of these peaks 
as a function of temperature. 

In summary, we investigated the plasmon modes of the 
Q1D electron systems on the surface of the liquid helium. 
We used a multisubband approach and treated the sys- 
tem as a classical nondegenerate gas obeying Maxwell- 
Boltzmann statistics. We found strong intersubband 
plasmon modes related to the first two subbands. The 
single-particle excitations in the system is responsible for 
strong Landau damping at frequencies w = noj n where 
gaps appear in the plasmon spectra. 
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FIG. 1. Spectral weight of the single-particle excitations 
S*p en {q,Lu) (solid curves) and S°* d (q,u) (dashed curves) for 
T = 0.4 K and N e = 10 4 cm -1 with different wavevectors q 
as indicated. 



[1] M. W. Cole, Rev. Modern Physics 46, 451 (1974); "Two- 
Dimensional Electron Systems in Helium and Other Sub- 
stracts ", Ed. E. Y. Andrei, Dordrecht: Kluwer (1997). 

[2] C. C. Grimes and G. Adams, Phys. Rev. Lett. 36, 145 
(1976). 

[3] Yu. Z. Kovdrya and V.A. Nikolaenko, Fizika Nizkikh 

Temperature 18, 1278 (1992). 
[4] J. Voit, Rep. Prog. Phys. 58, 977 (1994). 
[5] Yu. Z. Kovdrya, Low Temperature Physics 29, 77 (2003). 
[6] P. Glasson, V. Dotsenko, P. Fozooni, M.J. Lea, G. Papa- 

georgiou, S.E. Andresen, and A. Kristensen, Phys. Rev. 

Lett. 87, 176802 (2001). 
[7] S. S. Sokolov, G. Q. Hai, and N. Studart, Phys. Rev. B 

51, 5977 (1995). 
[8] S. S. Sokolov and N. Sudart, Phys. Rev. B 57, R704 

(1998); J. Phys. Condens. Matter 12, 9563 (2000); S. S. 

Sokolov and O. I. Kirichek, Low Temp. Phys. 20, 599 

(1994). 

[9] D.K. Lambert and P.L. Richards, Phys. Rev. B 23, 3282 
(1981). 

[10] Marcos R.S. Tavares and G.-Q. Hai, J. Phys.: Condens. 
Matter 13, 6421 (2001); G.-Q. Hai and Marcos R. S. 
Tavares, Phys. Rev. B 61, 1704 (2000). 



FIG. 2. The spectral weight (a) S cvcn (q, uj) and (b) 
S odd (q,u) for q = O.lfeo and T = 0.4 K. The electron den- 
sities are from iVe = 1 x 10 3 cm -1 (the lowest curve) to 
N e = 1.45 x 10 4 cm" 1 (the top curve) with a difference 
1.5 x 10 3 cm -1 . The insets indicate the position of the peaks 
as a function of N e . 

FIG. 3. The spectral weight (a) S cvcn (q,co) and (b) 
S° dd (q,u) for T = 0.4 K, N e = 10 4 cm" 1 , and different 
wavevectors corresponding to those in Fig. 1. The insets give 
the dispersion relation of the plasmon modes. The open-dots 
indicate the small peaks due to single-particle excitations 

FIG. 4. The spectral weight (a) S cvcn (q,uj) and (b) 
S° dd (q,u) for q = 0.1 k and N e = 10 4 cm" 1 . The tem- 
perature increases from T = 0.2 K (the lowest curve) to 2.0 
K (the top curve) with a step of 0.2 K. The insets show the 
position of the peaks. 
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